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A realistic model of a neutron star in minimal dilatonic gravity 
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We present a derivation of the basic equations and boundary conditions for relativistic static 
spherically symmetric stars (SSSS) in the model of minimal dilatonic gravity (MDG) which offers 
an alternative and simultaneous description of the effects of dark matter (DM) and dark energy 
(DE) using one dilaton field il?. The numerical results for a realistic equation of state (EOS) MPAl 
of neutron matter are presented for the first time. The three very different scales, the Compton 
length of the scalar field A<j, the star’s radius r*, and the finite radius of the MDG Universe ru, are a 
source of numerical difficulties. Owing to the introduction of a new dark scalar field ip = ln(l + ln il?), 
we have been able to study numerically an unprecedentedly large interval of A# and have discovered 
the existence of « 2.1 km for a neutron star with MPAl EOS. This is related to the bifurcation 
of the physical domain in the phase space of the system. Some novel physical consequences are 
discussed. 
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I. INTRODUCTION 

One of the most important lessons from the spectacular development of cosmology in the last fifteen years is the 
clear understanding that Einstein’s General Relativity (GR), as a model of gravity, and the Standard Particle Model 
(SPM), as a model of matter, are not enough to explain all the observed phenomena in Nature. 

There exist three possible ways for further development: a) To add some new content to the Universe, such as DM 
and DE; b) To change the theory of gravity^; c) The current observational data do not exclude a combination of a) 
and b). 

While the need for DM and DE has already been firmly established [l| , their nature and their small-scale distribution 
are still largely unknown. The only settled part concerning DM is its gravitational interaction. We have no evidence 
that DM has any other interaction but gravitational. We also have no idea what is the nature of DE. Within the 
framework of MDG [1^ DE can be thought of as responsible for two quite different physical processes with similar 
manifestations: the initial inflationary expansion and the present accelerating expansion of the Universe. 

The main general problem for the construction of f{R) theories still remains the absence of physical intuition when 
we are trying to specify the functio n f{R) [^. Nowadays, a series of additional requirements have been formulated 
aiming at cosmological applications |23l. l34|. However, in the literature one can find dozens of such functions. Several 
of them, for example BlilQIil, are thought of as valuable, bearing in mind cosmological applications. 

The situation in star physics is similar. The development so far, which has already lasted several decades, has not 
solved the problem of finding the real EOS of dense matter. At present, one can find several dozens of them dubbed 
realistic in the literature, see for example the very recent review (3 |. 

There have also been a series of attempts to use /(i?) models of gravity adapted to star physics^. For some constraints 
on f{R) for star models and specific numerical solutions, see No convincing final result has been reached. In 

the recent review one can find a description of the present state of affairs: While it is hard to construct Neutron 
Star (NS) equilibrium configurations in /(i?) gravity from a numerical point of view, there is no fundamental obstacle 
to their existence. NS configurations with realistic values of the physical parameters have never been constructed in 
viable f{R) models. 

The main goal of the present paper is to create a clear physical, analytical and numerical basis for the application 
of one of the simplest modification of GR, namely MDG, and to present for the first time a realistic model of NS in it. 

We hope that this consideration may help future developments of MDG models of the (almost) spherical objects 
at very different scales: laboratory scales, compact star scales, and at the scales of planets, white dwarfs, standard 
stars, star clusters, dwarf sphericals, galaxies, and clusters of galaxies. The use of the available information for similar 
physical phenomena at all reachable scales will give a much more definite justification of the model. A simultaneous 
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^ The two simplest modifications are the f{R) model, see for example O- Issll and MDG [36l - [S . These are similar, but not identical. Much 
more sophisticated models of modified gravity are also under investigation at present, see for example the references in [281 Isol . [43l . All 
of them have a lar ge n umber of additional parameters without clear physical justification and will be not considered here. 

^ See I 25 L l28l . [23.152l - f75l| and the references therein. 
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and coherent adequate description of these phenomena is a promising way to overcome the existing problems that are 
related to the absence of sufficient experimental and observational information about both matter and gravity. Such 
an ambitious program needs a well-developed theoretical, analytical and numerical basis. The best way to start is to 
probe well-known and simple realistic examples of SSSSjlike MDG-NS 

The MDG model was proposed and studied in (H, Is^, 113 • It describes a proper simple modification of GR 

based on the action A = Ag^^ -f Amatter- The action of the gravi-dilaton sector is 

= -^J d^x^\{^R - 2A17($)). (I.l) 

We take the Einstein constant k = 8ttGn/c^ ~ 1.8663 x cm g~^, Newton’s constant Gn ~ 6.6738 x 

10~^ and the cosmological constant A m 1.0876 x 10“®® cm“^. The dilaton field is $ G (0,oo). In 

general, this model is only locally equivalent to the f{R) model [13, and has a clear physical meaning: 

• The scalar field $ is introduced to take into account a variable gravitational factor G(<1>) = Gjv/d* = G]s[g{^) 
instead of the gravitational constant Gn and does not enter into Amatter, having no interaction with SPM 
matter. 

• The cosmological potential 17($) is introduced so as to have a variable cosmological factor A{7(d>) instead of the 
cosmological constant A. 

In GR with cosmological constant A, we have d) = 1, 5 (d)) = 1, and U{1) = 1. Due to its specific physical meaning, the 
field d> has quite unusual mathematical and physical properties and does not enter into the standard action Amatter ■ 
For astrophysical reasons, the cosmological potential !/($) must be a positive single valued function of d> G (0, 00 ). 
In [ 43 , there was introduced the class of withholding potentials, in order to confine dynamically the values of the 
dilaton d> in the physical domain, excluding antigravity, ghosts and tachyons. It is hard to formulate an analogous 
general requirement for f{R). Thus, we immediately see the main advantage of the MDG model, even when it is 
formally equivalent to some specific f{R) theory: We have a great deal of physical experience working with potentials 
like the cosmological one, both at the level of classical and quantum mechanics, or classical and quantum field theory. 
In units where Gn = c = 1, the field equations of MDG can be written in the form®: 

+ STrff = 0 , □$ + AV^ ($) = ^T. ( 1 . 2 ) 

The relation (d?) = | (d? (d>) — 2U (d>) j defines the dilatonic potential V (d>) (For the conventions used, see 

113). 

The main physical problem with all modifications of GR with one (or more) additional scalar field d> is the value 
of its mass m$. In the simplest cases of modified gravity, this is the only new parameter. This problem appeared for 
the first time as early as in (50| and is still remains unsolved. 

In Starobinsky (1980) presents an f{R) theory with one additional parameter^ /(R) = R + 6 m|,. Despite the 

fact that this model is still the most successful model of initial inflation [43, its potential R(d>) ^ ($ — 1 )^ allows 
antigravity and makes the model unacceptable in different physical situations. As longas one is working in a small 
enough vicinity of the vacuum state ^vac = 1, this shortcoming may be ignored at the classical level. In general, it 
will be not ignorable at the quantum level. Gomparing the scalaron mass with the cosmological data about the initial 
inflation, Starobinsky was able to find the value ^ 3 x 10~^ Mpiank [13 ■ This gives an extremely small Compton 
length A$ ^ 10“^^ cm of the scalar field. Such a model is indistinguishable from GR at the scales which we discussed 
above because of the Yukawa character of the corrections (~ ea;p(—r/A$)) to Newtonian gravity (^. [^. 1^. IsoL IsH . 

On the other hand, in the linear approximation, MDG reproduces the Yukawa tails to Newton’s law. The compar¬ 
ison with the laboratory and solar system experiments known in 1999 led to the estimate > 10“®eR/c^, which 
corresponds to A$ ^ 10“^ cm (^ . Later on, similar results were reproduced and refined many times in the framework 
of f(R) theories, see for example the recent papers [Hlal. According to the latest review of the experimental data 
1^ . the Compton length of the scalar field is A$ < 2.3 x 10“® cm. 

The simplest withholding dilaton potential in MDG can be written in the form^ [33I13|43- 

^ (‘J’+ - 2), (1.3) 


^ Here is the standard conserved energy—momentum tensor of the matter, Xa = X^ — ^XS^ is the traceless part of any tensor X^ in 
four dimensions, X = X^ is its trace, and the comma denotes differentiation with respect to 

The potential Ell has an unique minimum (the de Sitter vacuum) at ^vac = 1 and is the only withholding potential for which the 
corresponding Newtonian-like equation in a flat space time is solved by elliptic functions and the corresponding Schrodinger-like equation 
is solved by Heun’s functions. All other withholding potentials in these two cases will require the use of hyper-elliptic functions or of the 
not well-studied Fuchsian functions with more than four singularities. 
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and used in NS physics in Here p = A$-\/A is the dimensionless Compton length of the scalar field in 

cosmological units. According to the above estimates, this extremely small quantity lies in the physical interval 
p G (1 X 10“^®, 2.4 X 10“^^). Taking into account that the largest value was obtained only in the linear approximation, 
one can not be sure that much larger values of p are excluded by observations. Actually, to explain the current 
accelerating expansion of the Universe in the framework of the quintessence models, or the modified gravity models 
equivalent to them, one needs a very small mass of the scalar field ~ 2 x 10“^^ eV [t^. This gives A$ ^ 10^® cm 
and p ~ 10“^. The last value is still admissible in MDG (^1^. Hence, to check the MDG model at different physical 
scales, we need special techniques to deal with the interval p G 10“^). This is a very challenging task®. 

In MDG we can overcome the physical problem of the existence of different masses by introducing the potential 
U($) with many minima (^ . Then, around each of these minima, a Taylor series expansion will produce different 
effective masses of the scalar field. If so, the first step will be to study problems with simple potentials (11.31) for 
different values of the parameter p, and then to try potentials with many minima. One can hope that different effective 
values of will correspond to the above problems with different physical scales. 

In the present paper we study the problem at the scale of realistic NS. One of the goals is to recover the reasons for 
the numerical difficulties and to develop new methods to surmount them. 

We solve some of the physical problems by introducing the notions of cosmological energy density and pressure, 
and dilatonic energy density and pressure, see Eqs. (III. 21) and (III. 31) . as well as novel equations of state for them: 
Cosmological EOS (CEOS) and Dilatonic EOS (DEOS), which have to be used in MDG together with the familiar 
Matter EOS (MEOS) An essential role is played by the finiteness of the dilaton pressure at the center of a 

star, see Eqs. (III. 41) . Here we show for the first time in detail how these notions and relations arise from MDG (See 
Appendix lAI). 

A star is a matter excitation above the proper vacuum state of the theory. Since in MDG the physical vacuum 
is the de Sitter vacuum we define the surrounding mass m(r) of the star using the relations —Ijgrr = 

1 — 2m{r)jr — |Ar^ (4H448| . Note that, in those articles on stars in f{R) models known to the author, the mass m{r) 
is often defined using the incorrect assumption that the vacuum state corresponds to a fiat spacetime asymptotic, i.e., 
without the A term. 

Outside the star, the dilaton creates a “dilaton sphere,” or, for short, a “disphere” In the non spherically- 

symmetric case, one can speak of a “dark domain” outside the matter. The disphere exponentially dilutes up to the 
radius of the Universe ry, defined by the cosmological horizon of the MDG model where the de Sitter vacuum is 
reached. Here, we calculate, for the first time, the mass of the disphere of NS with a realistic MEOS, and the total 
mass rritotai of the whole complex: NS plus disphere. Some similar results and a not very precise use of the terminology 
of articles can be found in devoted to quark stars with MIT-bag-model MEOS. 


II. BASIC EQUATIONS AND BOUNDARY CONDITIONS FOR SSSS IN MDG 


The spacetime interval for SSSS is ds^ = — r^dVi^, where r is the luminosity distance to the 

center of symmetry and dVi^ describes the space-interval on the unit sphere®. After some algebra (See Appendix lAl) 
one obtains the following basic results for SSSS: 

In the inner domain r G [0,r*], the SSSS structure is determined by the fourth order generalization of the TOV 
equations: 

dp p + e m + Anr^pef //$ 
dr 

dp^ _ p^ (^ _ 2.c( , q ,^A 2 

dr rA 


A 2 

— = Airr ee///$. 


dd) 

dr 


= -Airr'^p^/A, 


A — 27rr^p^j^ 

^3r — 7m — + 47rr®ee///$^ — 


(H.la) 

(H.lb) 


The four unknown functions are m{r), $(r), p^ (r), and p{r). In Eqs. (III. II) A(r) = r — 2m(r) — ^Ar®, Ce// = e-l-e^ -|-e^, 
Pef f = P+Pa AP* ■ addition to the standard MEOS e = e{p), we obtain two novel equation: CEOS: 
and DEOS: e, = p - ie + ($) + ^ 

The cosmological energy density and the cosmological pressure are defined as follows: 

A(c/(*)-<t). („.2) 




^ One way to get around this problem is to use some sophisticated additional mechanisms like the chameleon, K-mouflage, or Vainshtein 
mechanisms, see for example |77L ItSI and the references therein. Then the mass depends on the environment in a somewhat ad hoc 
and artificial way. 

® The luminosity radius r is an invariant defined by the relation A = ; A is the area of a sphere around the center of symmetry. 
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The dilatonic energy density and the dilatonic pressure measure the changes of the gravitational factor. By defini¬ 
tion:^ 




V~9r 


47rr dl 


A 

Airr^ 




Placing the physical SSSS-center at Tc = 0, we obtain the boundary conditions (See Appendix iBll: 


(11.3) 


m(0)=mc = 0, 4>(0) = 4>c, p{0) = Pc, P* (0) = ^ (4>c). (II.4) 

Requiring rric = 0 ensures the finiteness of the pressure Pc simultaneously for the Newton-, GR- and MDG-SSSS. 
The condition on p$c (= ~|e<i>c) ensures its finiteness, being a specific relation for the values of the MDG-center: 
Rb(p$c,Pc,^’c) = 0. 

The SSSS-edge is defined by the condition p* = p{r*;pc, 4>c) = 0. Then 

TO* = $* = $(r*;pc, $c), pI = P^ir*-,pc,<^c)- (11.5) 

The luminosity radius of a compact star with physically realistic MEGS may vary: r* ft! 8 14 km. 

To obtain a complete description of the spacetime geometry inside SSSS, one must add Eq. (IA.3I) . which splits out 
from Eqs. (|II.1I) . as well as the corresponding boundary conditions for i^(r). 

Outside the star, p = 0 and e = 0, and we have a disphere (45l - l48l |. Its structure is described by the shortened 
system (nnj, where the second of Eqs. (III.la|) is omitted. Eor the exterior domain r G we use Eqs. (III. 51) 

as left boundary conditions. The right boundary is defined by the MDG cosmological horizon Tjj-. A{r^-,pc, 4>c) = 0, 
where the de Sitter vacuum is reached: ^(r^;pc,^c) = 1 and gu = 0, prr = OO) gttgrr = 1. As a result, we obtain a 
new relation for the values of the MDG SSSS center, Fa{pc, ^c) = 0 which depends also on the Compton length A$. 
Besides, the 4d and 3d scalar curvatures are = 4A — 1/r^ and = 4A. A schematic picture of MDG Universe 
with only one SSSS is shown in Fig. [T] 
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FIG. 1: A schematic picture of MDG Universe with one SSSS 

The two MEOS dependencies ‘kc) = 0, Ea(pc, ‘kc) = 0 show that for a given MEOS in MDG, as well as 

in Newtonian gravity and GR, we have a one-parameter family of SSSSs. 

Now it becomes clear that a basic difficulty in the numerical investigation of NS in MDG is the presence of three 
very different scales: A$ G cm, 10^^ cm), r* ^ 10® cm, and rjj ^ 10^® cm. The justification of the values of A<i> 

for different astrophysical objects becomes the most important physical problem for modified gravity. 

In the MDG Universe, the maximum mass of any matter object is Mmax = lO^^M©. This value is physically 

safe as it is about six orders of magnitude greater than the mass of the most massive objects in the Universe [79|. 
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Here we use the true geometrical radial distance l(r) = 2m(r) a monotonically increases from 0 to some finite value 

V 1 F "3 ^ 

> '^U (the size of the Universe), since ^ > 1; l{r) r, when r — >■ 0; and — >■ oo, when r —>■ r^/. 
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III. A MODEL OF A NEUTRON STAR WITH MEOS MPAl 

We use the well-known realistic MEOS MPAl (sol - l^ . Its analytic version [s^ allows a treatment of NSs with 
central densities up to 10^® g/cm^. The matter density decreases to Ri 6.5 g/cm^ on the surface of the NS. 

To have successful computations, we were forced to implement high-precision computer arithmetic with 64 digits 
and to replace the dilatonic field $ with a novel field variable 

p = ln(l -p ln$) $ = exp(exp((p) — 1). (IH-l) 

We call the field the dark scalar. Note that the mass of the dark scalar precisely equals the mass of the dilaton m$. 

The double logarithmic substitution (llll.f ll stretches the physical domain of the scalar field and makes possible 
numerical calculations in the maximal interval for the parameter p allowed by the given MEOS. Inside the star we use 
standard logarithmic variables ^ = log(p) and ( = log(p). Outside the star we also use a proper logarithmic variable 
instead of the luminosity radius r. The numerical results were obtained by an appropriate version of the shooting 
method. 

The key step in the calculations is to obtain the relation Ea((/3c, Cc! = 0 shown in Fig. [5] for different values of 
the Compton length A$. The boundary conditions at the center of the star defined by these curves yield the MDG 
iTT-totai — T* relations shown in Fig. |31 The dashed black line describes the corresponding OR m* — r* relation. As 
it should be, in the limit A$ —>■ 0 the MDG rritotai — r* curves tend to the GR one. Fig. U] shows the corresponding 
dependencies of the compactness C = m* jr* of an MDG NS on the central matter density. 



FIG. 2: The specific MDG curves Fa{Pc,^c', As) = 0 for NS with MPAl MEOS and different fixed A# £ (2.4 km, 10“^ km). The 
curves, from top to bottom, correspond to decreasing values of A*. 



r* km 


FIG. 3: The MDG mass-radius relation for a NS with MPAl MEOS for different hxed 
A# £ (2.4,10^) km. The black dashed curve presents the GR mass-radius relation 

As seen in Figs. [2H11 the influence of the dark scalar on the interior structure of the NS is significant. The decrease 
of the value of A$ leads to a narrowing of the domains of the corresponding variables. When A$ approaches the critical 
value the curves shrink to a point. This is an indication of the existence of a bifurcation point of the physical 
part of the phase space of the system. For an NS with MPAl MEOS we obtained numerically ~ 2.1 km. 

In Figs. [SHU one can see more consequences of the presence of the dark scalar in NS physics. 

After all, as seen in Fig. 0 the MGD NSs have qualitatively the same stability properties as the GR ones. Clearly, 
the numbers depend on the dark scalar mass and may vary to some extent. 
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FIG. 4: The compactness of an MDG NS with MPAl MEGS for different fixed 
A$ £ (2.4 km, 10“^ km) That for GR is indicated by the black dashed curve. 



r km 



FIG. 5: Examples of dependence of gravity intensity on r. The dark scalar, having no direct interaction with the matter of 
SPM, influences the structure of the NS by changing significantly the intensity of gravity inside it, and in its vicinity. 




FIG. 6: Left: The dependence of the surrounding mass m(r) on the luminosity radius r. Right: mtotai versus m* 



C 


FIG. 7: The mtotai - Cc dependence of MDG NS with MPAl MEOS 


IV. DISCUSSION 

The present paper is devoted to MDG as the simplest possible solution of the current problems with DM and DE. 
Here we considered for the first time a realistic MDG model of a neutron star (NS) with MPAl MEOS reaching the 
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following basic results: 

The derivation of the MDG generalizations of the TOV equations and the corresponding boundary conditions at 
the star’s center, at its edge, and at the cosmological horizon of the de Sitter like MDG Universe. A clear physical 
understanding of the SSSS structure is reached using notions of cosmological energy density and pressure, dilatonic 
energy density and pressure, and the corresponding equations of state: CEOS, DEOS, and MEOS. 

The maximal mass M^ax = lO^^M© of any matter object in an MDG Universe is consistent with observa¬ 

tions. 

The existence of three very different MDG scales, the Compton length of the scalar field A$ G 10^®) cm, 

the star’s radius r* ~ 10® cm, and the finite radius of the MDG Universe ru ^ 10^® cm, is a source of numerical 
difficulties. Owing to the introduction of a new dark scalar field ip = ln(l -|-ln$) we were able to numerically study for 
the first time an unprecedentedly large (four orders of magnitude) interval of A$ for neutron stars with MPAl MEOS 
and to discover the existence of A^’'** « 2.1 km, related to the bifurcation of the physical domain in the phase space of 
the system. This value corresponds to a critical mass r; 5 x 10“^^ eV j<? and depends on the MEOS of the NS. 

The kind of narrowing, typical for a bifurcation point, of the domains of the corresponding variables (see Figs, 
may explain the astrophysical observations which do not show the existence of NSs along the whole theoretical mass- 
radius curves, but only on a narrow part of them (ssl - l^ . This unexpected result of ours may serve as a new criterion 
for the choice of a realistic MEOS, coherently and simultaneously with the determination of the dark scalar mass to$. 
For this purpose, one needs a more profound application of the mathematical theory of bifurcations and we intend to 
publish the corresponding results separately. 

The gravitational force in the interior of an MDG star is smaller than in GR and may vary in different ways, see 
Fig. [51 This result refutes the wide-spread opinion that the A term in the gravitational action is inessential at the 
stellar scale. Indeed, in the physical de Sitter vacuum, we have $ = 1 and U = 1 and the extremely small observed 
value of A makes negligible the A-terms at the scale of the solar system and in a large enough vicinity outside a star. 
However, inside the star, we have no vacuum state and the dilaton deviates from its vacuum value. As a result, the 
cosmological potential changes its value very significantly inside the star and compensates for the small value of the 
cosmological constant. 

The above statements and the previous work on MDG opens a novel possibility: To look simultaneously and 
coherently for a realistic MEOS of different physical objects and for a realistic withholding cosmological potential, 
which together are able to describe the variety of phenomena at very different physical scales. 

It would be very interesting to work out models of moving and rotating stars in MDG. In this case, one can expect 
not only an asymmetric stellar configuration and dark domains, but also the appearance of different centers of the 
star and its dark domain, or even detachment of par ts of the dark domain. At the galactic and galactic cluster scales, 
such phenomena have already been observed (88l - l^l | and may allow an MDG explanation. 
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Appendix A: The field equations for static spherically symmetric MDG 


Taking into account the expression for the scalar spacetime curvature R=—e ^ 
from the basic Eqs. den one obtains the system of ordinary differential equations 
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+ {I-I (3-' + A')) (-' - A') - " (-' + A')) $ + ( 24 . {e + p)-^<4>'^ 

-3$" + 0 + i (i.' + 3A')) ‘^'+(^^-^''-1^' (^' - A') + ^ W + A')) - (s^ (^+P) + ^^) 

- A')) ® (e +P) - 

$" 2 ^^' ~ A^^ + (Stt ^- — e^ 
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0, (A.la) 
0, (A.lb) 
0, (A.lc) 
0. (A.ld) 


The first three of Eqs. (lA.ll ) are (*). (^). and (^) proiections of the first of equations in Eqs. (iLil, and Eq. (|A.ld|) 
follows from the second equation of (11.21) . 

Because the first equation of Eqs. m is traceless, Eqs. (IA.lal) - (|A.lcl) are not independent and one can omit Eq. 
(lA.lcE Then, using the relation R = 217 $ [d^, the explicit form of the scalar curvature R, and the definition of y$, 
one obtains from Eqs. (lA.lal) and (lA.lbl) 


Let us introduce the mass function m(r) obeying the relations = A/r, A(r) = r — 2m{r) — The next step is 
to define the dilaton pressure and dilaton energy density according to relations (III. 31) . As a result, Eq. (|A.2al) acquires 
the form (Ill.laE Using Eqs. (in.3p we obtain also the second equation of (III.lbE 

Using the definitions of Peff and pA, from Eq. (IA.2bl) we obtain an equation for the function i/{r) that is valid both 
inside and outside SSSS: 

7 = (A,3) 

r A — 27rr^p$ /$ 

Outside the star, the function i/(r) is determined by Eq. (IA.3I) under the additional conditions p = 0, e = 0. 

Inside the star we use the standard approach: in MDG, as well as in GR, the equation of motion of matter is Vj/T^ = 
0. We assume that SSSS is filled by a standard ideal fluid with energy-momentum tensor = {e+p)u^u'^ —p5'^^ at rest, 
i.e., when U/j, = (5°. Then, in thermodynamical equilibrium, inside the star we have the usual relation p' = —{e+p)i''/2. 
Inserting Eq. (|A.3I) into it, we obtain the second of Eqs. (III. lap , thus arriving at the extension (III.ip of the TOV 
equations in the case of MDG SSSS. 

Finally, in terms of these newly introduced quantities, Eq. (lA.ldp transforms into DEOS. 


Appendix B: Conditions at the center of the star and the behaviour of the solutions 

As in Newtonian gravity, and in the GR gravity of SSSS, we assume that all physical quantities are finite at the 
center Tc = 0 of the MDG star. The withholding property of the cosmological potential [/(<!>) dynamically ensures 
that 0 < <&£ < 00 [i^l and, as a result of their definitions, Uc, Cc, V% , and p\ are finite. As usual, we assume 
0 < Ec < 00 and 0 < Pc < 00 for the stellar matter. In addition, we require |p|| < 00 and |e|,| < 00 . Then Ip^ffl < 00 
and < 00 . 

Now, taking the limit r0, we obtain from the first of Eqs. (Ill.lal) m' ^ => 77i(r) ^ TOc +|7rr^egj^/<I>c- 

If the integration constant rric ^ 0, then the second of Eqs. (Ill.lal) yields p(r) ^ ^ (Pc + Cc) ln(r/ro) + const. Since 
Pc + Ec > 0, then p(r) —>■ —00 for r —>■ 0. This is physically unacceptable, just as in Newtonian gravity and in GR. 
Thus we see that one must suppose that TOc = 0. Then 

m(r) ^ + 0{r^), p{r) + e,) + 0(r3). (B.l) 

From the second of Eqs. (Ill.lbl) we obtain p^ ^ —yP| — P*(^) — (3p| + 2e|) ln(r/ro) + const. Obviously, 

the only physical solution with |p|, | < 00 is the one with 

P$ = ~|e$ = gTc — 77 = Cc — 3 pc. 


(B.2) 
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Using this result, from Eqs. (Ill.lbl) we obtain 


$(r) - + fTrr^el + 0 {r^) = 



m{r) 


■ff r 



(B.3) 
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